Stabilization of the p-wave superfluid state in an optical lattice 
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It is hard to stabilize the p-wave superfluid state of cold atomic gas in free space due to inelas- 
tic collisional losses. We consider the p-wave Feshbach resonance in an optical lattice, and show 
that it is possible to have a stable p-wave superfluid state where the multi-atom collisional loss is 
suppressed through the quantum Zeno effect. We derive the effective Hamiltonian for this system, 
and calculate its phase diagram in a one-dimensional optical lattice. The results show rich phase 
transitions between the p-wave superfluid state and different types of insulator states induced either 
by interaction or by dissipation. 
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The observation of the s-wave superfluid state in a 
Fermionic atomic gas represents a remarkable break- 
through in the study of many-body physics with ultra- 
cold atoms l]. The Feshbach resonance plays an impor- 
tant role in those experiments, enhancing the interatomic 
interaction so that the superfluid phase can be entered 
at a temperature that is experimentally achievable. It 
is of great interest to realize the superfluid state with 
other pairing symmetries as well. The p-wave superfluid 
state is the next candidate for observation, and has at- 
tracted the interest of many experimental and theoretical 
groups [3, 0] ■ The p-wave Feshbach resonance has been 
recently observed in experiments and can push the 
single-component Fermi gas to the strongly interacting 
region and open a door towards observation of the p- 
wave superfluid state in this system. However, compared 
with the s-wave Feshbach resonance, a key difficulty with 
the p-wave resonance is that the inelastic collision loss in 
this system is typically large @, 0; Q : which forbids ther- 
malization of the gas within the system lifetime [sj . 

In this paper, we dicuss how a dissipation-induced 
blockade mechanism can stabilize the p-wave superfluid 
state in an optical lattice in the strongly interacting re- 
gion. Such a dissipation-induced blockade has been re- 
jorted in recent experiments to realize the Tonks gas 
6 1 or for simulation of effective three-body interactions 
7| with cold bosonic atoms. We apply this mechanism 
to stabilize the single-component Fermion system in an 
optical lattice in the presence of the p-wave Feshbach res- 
onance. The p-wave Feshbach resonance has been con- 
sidered recently on single sites in a very deep optical lat- 
tice 0, Here, instead, we focus on the many-body 
physics by deriving an effective Hamiltonian for this sys- 
tem, taking into account the atomic hopping in the re- 
duced Hilbert space caused by the dissipation-induced 
blockade. This effective Hamiltonian provides a starting 
point to understand the quantum phases. We compute 
the phase diagram of the system explicitly with well- 
controlled numerical methods for an anisotropic lattice 



where the atom tunnelling is dominantly along one di- 
mension. The results show rich phase transitions between 
the p-wave superfluid state, a dissipation-induced insula- 
tor state, the Mott insulator state, and different kinds of 
metallic states. Although these results are obtained from 
one-dimensional calculations, we expect these phases to 
correspond also to similar phases in higher dimensions. 

We consider a single-component Fermi gas near a p- 
wave Feshbach resonance. If this strongly interacting gas 
is loaded into an optical lattice, many different Bloch 
bands can be populated, in particular when the resonance 
is broad (as it is the case for some recent experiments 0). 
However, we can derive an effective single-band model for 
this system that is independent of the interaction details. 
Following a strategy similar to the s-wave Feshbach res- 
onance case [9], we first analyze the local Hilbert space 
structure on a single lattice site. When we have zero or 
one atom on the site i, the states are simply denoted by 
|0)- and |a)- — a\ |0)-, respectively. For the case of two 
atoms on a single site, the exact two-body physics has 
recently been calculated [Sj], and there are several two- 
body energy levels separated by an energy difference of 
the order of the band gap. If we assume that the system 
temperature is significantly below the band gap, only the 
lowest two-body state is relevant. We refer to this state as 
a dressed molecule level, and denote it by \b)^ = b\ |0)j. 
Note that the wave function of \b)^ in general includes 
contributions from many of the original atomic orbitals 
Q. It has a p-wave symmetry in space and is antisym- 
metric under exchange of the two atoms in this dressed 
molecule. 

If more that two atoms coming to a single site, different 
from the s-wave case, the state will not be stable due to 
big three-particle inelastic collision loss At first sight, 
this seems to mean that the system will become unsta- 
ble. However, in a lattice, there is a subtle dissipation in- 
duced blockade mechanism Q which forbids population 
of the unstable three-particle state and thus stabilizes 
the whole system. The basic idea is illustrated in Fig. 1. 
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FIG. 1: (Color Online) Illustration of the dissipation-induced 
blockade for multiple occupation (more than two) of a single 
lattice site. Similar to a two level transition with a detuning 
Aa and a decay rate 7 for the target level, the effective hop- 
ping in this case is suppressed by a factor t^/ (7^ -I- A|), where 
t is the atomic tunneling rate. 



The three-particle state |3) has a large bandwidth char- 
acterized by its inelastic collision rate 7 and an energy 
shift characterized by the on-site atom-dressed-molecule 
interaction energy A3. If a single atom tunnel through a 
barrier with a hopping rate t to form this state |3), the 
probability of getting |3) is given by / (7^ + A3) (sim- 
ilar to a two-level transition with a detuning A3 and a 
bandwidth 7). So the net colhsional loss of the system is 
bounded by 7off = -ft^ / (7^ + A|) < / (2A3) no mat- 
ter how large the inelastic collision rate 7 is. Near the 
Feshbach resonance, the atom-dressed-molecule interac- 



tion energy A3 is comparable with the lattice band gap 
(thus much larger than t) [llj , the net colhsional loss 7eff 
is therefore small compared with the atomic hopping re- 
quired to thermalize the system [ld|. The reduction of 
population in the noisy state |3) is called the dissipation- 
induced blockade (or interpreted as the quantum Zeno 
effect [1]; the blockade is actually induced by both dissi- 
pation and interaction when A3 and 7 are comparable). 
Due to this mechanism one can achieve many-body ther- 
mal equilibrium in a lattice even if there is big inelastic 
collision loss. 

Due to the dissipation-induced blockade discussed 
above, on each site we have only three relevant levels: 
0)j , |a)^ , and as the low energy configurations. The 
energy difference between the sates \a)^ and \b)^ can be 
tuned with the external magnetic field via the Feshbach 
resonance. We then take into account the atomic hop- 
ping, which changes the level configurations of the neigh- 
boring sites under the constraint that the atom number 
be conserved. Thus, there are only three possible pro- 
cesses, as illustrated in Fig. 2. Corresponding to these 
configuration changes, the general Hamiltonian for this 
lattice system then takes the form 



H = J2 - ^M<^^ + 2&IfeO] ~Y.^ [t^ataj + t2{b+ - b+)a,aj + t^b+hja+a, + H.c] P, 



(1) 



where /i is the chemical potential and A is the energy de- 
tuning of the dressed molecule controlled with the mag- 
netic field. The value of A characterizes the on-site 
atomic interaction magnitude. The hopping rates for the 
three processes illustrated in Fig. 2 are different, in gen- 
eral, due to the contributions from different bands (be- 
cause the dressed molecule is a composite particle with 
population in multiple bands). The hopping takes place 
in the low energy Hilbert space specified by the projector 
P = (|0)j (0| + |a) ■ (a| + \b),-^ (&|), and the summation 

i 

in Eq. (1) is over all neighboring sites 

The Hamiltonian (1), although much simplified com- 
pared with multi-band models, is still complicated, and in 
general it does not allow exact solutions. To understand 
some basic physical properties of the system, we limit 
ourselves in the following to the one-dimensional lattice 
where the atomic hopping along the other two dimensions 
are turned off with a high lattice barrier. In this case, 
we can solve this model with well controlled numerical 
simulations, and the results show rich phase transitions 
between the p-wave superfluid state and different kinds 
of insulator and metallic phases. We expect most of the 
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FIG. 2: (Color Online) Illustration of different atomic hop- 
ping processes over the two neighboring sites described by the 
Hamiltonian (1), where each site has only three possible level 
configurations. 



phases found in the one-dimensional case have counter- 
parts in higher dimensions. In particular, the p-wave su- 
perfluid state characterized by a quasi-long-range pairing 
order with diverging pair susceptibility in one dimension 
can be easily stabilized to a true long-range order if we 
allow weak tunneling between the one-dimensional tubes 
12|. 

In the numerical simulation, we use the iTEBD algo- 
rithm, a recently developed method related to density 
matrix renormalization group techniques [l3 , 14 1 , which 
allows direct calculation of the physical properties in the 
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FIG. 3: (Color Online) The dressed molecule number nt (the 
double occupation probability) and the atom number ria (the 
single occupation probability) shown as the function oi A/t 
(with a fixed fi/t = —0.5) in Fig. (a) and fi/t (with a fixed 
A/t — —1.5) in Fig. (b). The non-analyticities of these curves 
signal a number of quantum phase transitions from the dis- 
sipation induced insulator state (DII), to a p-wave superfluid 
state (PS), to a Mott insulator state (MI), to a normal mixture 
state (NM), and finally to normal Fermi gas phase (NFG). 



thermodynamic limit. The algorithm has been shown to 
work with high precision compared to known results for 
the Hubbard model [ll] . For simplicity, we take the hop- 
ping rates ti = t2 — — t, and use t as the energy 
unit. Then we have effectively only two parameters, A 
and fji (in units of t) in the Hamiltonian (1). To figure 
out the complete phase diagram with respect to these 
two parameters, we calculate d{H) /OA and d {H) /dfi 
as functions of A or fi for the ground state of H, and use 
the characteristics of these curves to identify the phase 
transition points. 

In Fig. 3, we show nf, — ib\bi^ and ria = i^o^o^i^ 
as functions of A and ji. One can clearly see several 
quantum phase transitions from this figure. First, with 
a fixed chemical potential /i = — 0.5t, one has rlf, = 1 
and ria — with a large negative detuning A (corre- 
sponding to strong attractive atomic interaction) . In this 
case, each site is doubly occupied with two atoms. More 
than two atoms cannot move to the same site because 
of the dissipation-induced blockade. So this is an insula- 
tor phase stabilized by the dissipation. As one increases 
A with A > — 2.15f, the number of atoms on each site 
begins to fluctuate. If one looks at the pair correlation 



(blbjj, it shows quasi-long-range behavior with slow al- 
gebraic decay. In Fig. 4(a), we show this correlation 
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FIG. 4: (Color Online) The correlation functions shown for 
the p-wave superfluid phase ((a) and (b)), the Mott insulator 
phase ((c) and (d)), and the normal mixture phase ((e) and 
(f)). Figures (a), (c) and (e) show in fc-space the Fourier trans- 
form Pfc of the pair correlation (blbj^^. Figures (b), (d) and 

(f) show the charge density correlations of {nbiUbj), {nbiUaj) 
and {uairiaj) in real space. We take n/t = —0.5 for all the 
figures. 



in the fc-space, defined as Pk — (l/N) 



N 

E 

p=0 



The correlation Pk is peaked sharply at k — 0. This 
corresponds to the p-wave superfluid phase. The p-wave 
character is inherited from the p-wave symmetry of the 
dressed molecule in space bi (or the atomic pair on the 
same site) . The p-wave nature of the pairing is also man- 
ifested in the atomic pair wavefunction at different sites 
(fliaj), which is obviously antisymmetric under exchange 
of the sites. In the one-dimensional case, the p-wave su- 
perfluid state is characterized by a quasi-condensate of 
the atomic pairs with a diverging pairing susceptibility. 
If one allows weak coupling between the one-dimensional 
tubes in the optical lattice, the p-wave quasi-condensate 
can easily be stabilized into a real condensate with a true 
long range pairing order. 

If one further increases A in Fig. 3(a), one enters 
a phase where the total particle number per site n — 
2fib + ria is fixed at 1 (although the double occupation 

probability (blbi^ varies with A). This is a Mott insula- 
tor state with a finite gap to charge excitations. This can 
be seen more clearly in Fig. 3(b), where we fix the detun- 
ing A, and show and Ha as functions of the chemical 
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FIG. 5; (Color Online) The complete phase diagram of the 
system versus two parameters ^/t and A/i. The five different 
phases are marked with the same notation as in Fig. 3. The 
black dashed lines a and b correspond to the parameters taken 
in Fig. 3(a) and Fig. 3(b). 

potential ^. For this phase, the number density does 
not vary with so the system is incompressible as one 
expects for a Mott insulator phase. The correlation func- 
tions for this phase is shown in Fig. 4, where both the 
pair correlation and the charge density wave correlations 
are of short range. As one further moves to the right side 
in Fig. 3(a), there are two other phases: the normal mix- 
ture (NM) and the normal free gas (NFG). Both of these 
two phases are of a metallic nature with a finite com- 
pressibility (see Fig. 3(b)). The difference is that in the 
normal mixture phase, some sites are doubly occupied 
(with a finite (b\bi^). Several kinds of correlation func- 
tions for the normal mixture phase are shown in Fig. 4, 
and all of them decay rapidly with distance. In the NFG 
phase, the double occupation probability ib\bi^ reduces 

to zero, and one has a conventional single component free 
Fermion gas. 

We have calculated the phase transition points for all 
A and /x, and the result is summarized in Fig. 5 to give 
a complete two-parameter phase diagram. The five dif- 
ferent phases are marked there. Under real experimental 
conditions, there is typically a weak global harmonic trap 
in addition to the optical lattice potential. As one moves 
from the trap center to the edge, the effective chemical 
potential fx gradually decreases under the local density 
approximation. So with a fixed interaction parameter A, 
a line in the phase diagram of Fig. 5 across different /i 
gives the phase separation pattern of the Fermi gas in a 
harmonic trap. With a large positive A, one has a Mott 
insulator state in the middle, surrounded by a normal 
Fermi gas. The phase transitions are most rich for small 
I A I , where one can cross all the five different phases from 
the trap center to the edge. For large negative A, the 
region of the p-wave superfluid state increases, but the 
Mott insulator and the normal mixture states eventually 
disappear when A < — 3^. As the p-wave superfluid state 



has a large stability region in the phase diagram, such a 
phase can be prepared experimentally by adiabatically 
ramping the Hamiltonian parameters following certain 
trajectories that suppress the three-particle occupation 
0. 

In summary, we suggest that the p-wave superfluid 
state near a Feshbach resonance can be stabilized in 
an optical lattice through a dissipation-induced block- 
ade mechanism. We have derived the effective Hamilto- 
nian to describe this strongly interacting system, taking 
into account the restriction of the Hilbert space due to 
this blockade mechanism. We solve the Hamiltonian in 
the anisotropic one-dimensional lattice through exact nu- 
merical calculations, and the result suggests rich phase 
transitions between the p-wave superfluid state and sev- 
eral kinds of insulator or metallic phases. 
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